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Abstract: We study the power asymmetry between even and odd multipoles in the multipolar 
expansion of CMB temperature data from WMAP, recently reported in the literature. We intro- 
duce an alternate statistic which probes this effect more sensitively. We find that the data is highly 
anomalous and consistently outside 2a significance level in the whole multipole range I = [2, 101]. 
We examine the possibility that this asymmetry may be caused by the foreground cleaning proce- 
dure or by residual foregrounds. By direct simulations we rule out this possibility. We also examine 
several possible sub-dominant foregrounds, which might lead to such an asymmetry. However in 
all cases we are unable to explain the signal seen in data. We next examine cleaned maps, using 
procedures other than the one followed by the WMAP Science team. Specifically we analysed 
the maps cleaned by the IPSE procedure, Needlets and the harmonic ILC procedure. In all these 
cases we also find a statistically significant signal of power asymmetry if the power spectrum is 
estimated from masked sky. However the significance level is found to be not as high as that in the 
case of WMAP best fit power spectrum. Finally, we test for the contribution of low-Z multipoles 
to the observed power asymmetry. We find that if we eliminate the first six multipoles, I = [2, 7], 
the significance falls below 2a CL. Hence we find that the signal gets dominant contribution from 
low-? modes. 



1 Introduction 

The primary aim of WMAP satellite has been to measure full sky CMBR temperature anisotropies with 
great precision (Bennett et al., 2003a). The primary quantity of interest from these full sky maps is the 
temperature power spectrum, which is used to constrain various cosmological parameters (Hinshaw et al., 
2003; Spergel et al., 2003; Larson et al., 2011; Komatsu et al., 2011). It also gave TE cross power spectrum 
and E-mode power to a good precision (Kogut et al., 2003; Page et al., 2003; Larson et al., 2011; Komatsu 
et al., 2011). Since the release of WMAP 1st year data many large scale anomalies were reported in the 
data (Bennett et al., 2003b, Efstathiou, 2003; de Oliveira-Costa et al., 2004; Eriksen et al., 2004a,b; Ralston 
& Jain, 2004; Copi & Huterer, 2004; Schwarz et al., 2004; Land & Magueijo, 2005a; de Oliveira-Costa & 
Tegmark, 2006; Copi et al., 2007; Samal et al., 2008; Copi et al, 2010; Bennett et al., 2011). Recently, in 
(Kim & Naselsky, 2010), (see also de Oliveira-Costa et al., 2004; Land & Magueijo, 2005b; Gurzadyan et al., 
2007; Gruppuso et al., 2011; Hansen et al., 2011; Ben-David, Kovetz & Itzhaki, 2011), it was found that the 
CMB power in odd multipoles is anomalously more than that in the even multipoles. Some possible sources 
of such odd modulations to CMB, like signals from inflationary era and solar system physics, were discussed 
in (Groeneboom et al., 2010; Koivisto & Mota, 2011; Maris et al., 2011). This power asymmetry between 
even and odd multipoles is also addressed as Parity asymmetry. 

The CMBR temperature fluctuations on a sphere are usually expanded in terms of spherical harmonics, 
-Mm j as 
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where a; m 's are the multipolar expansion co-efficients. The power spectrum of CMB is defined as 
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We denote Q = 1(1 + \)Ci/2n. Kim & Naselsky (2010) defined the following two quantities 



1=2 



and 

2=2 

Here P + and P~ represent the sum of power in the even and odd multipoles respectively. The parity 
asymmetry statistic was defined as P + / P~ , which we refer to as "-P(Z m aa;)" for convenience. By comparing 
the data (Larson et al., 2011) with "pure" realizations of CMB signal generated based on KCDM model, 
they estimated the p— value to be 0.002. When this statistic was applied on full-sky Ci recovered from masked 
maps, using KQ85 mask from WMAP's seven year data release, the minimum probability or value was 
found to be p — 0.003. The l ma x corresponding to this p— values was found to be 22 from simulations, which 
is an a posteriori choice. By accounting for this posterior choice of l m ax = 22, the probability was estimated 
to have lowered to 0.02. 

In the present paper we study this parity asymmetry in considerable detail. We first consider an alternate 
statistic to test for parity asymmetry. This statistic is found to be more sensitive than the one considered 
in (Kim & Naselsky, 2010). We next investigate whether this asymmetry might arise due to foregrounds, 
which are not symmetric under parity. Most of the foreground contamination, however, gets removed in 
the process of extracting the primordial power spectrum. Nevertheless the residual foregrounds might be 
sufficiently large to cause parity asymmetry. Hence we use simulated foreground cleaned maps to test the 
significance of parity asymmetry in the WMAP CMBR data. We utilize both the ILC and IPSE cleaning 
procedure for this purpose. The simulated maps are generated using random realizations of CMBR and 
pre-launch Planck Sky Model for foregrounds. We also allow for the possibility of some unknown foreground 
components. Another interesting anomaly found in the CMB data is the ecliptic dipolar modulation of the 
CMB power, discovered in (Eriksen et al., 2004a). This signal is also parity asymmetric and hence one may 
suspect that there might be a relationship between this and the signal discovered in (Kim & Naselsky, 2010). 
We study the possibility that a dipole modulation of temperature anisotropy might lead to the observed 
parity asymmetry. Furthermore we examine whether foreground cleaned maps obtained using alternate 
procedures such IPSE, Needlet ILC, etc., also show the parity asymmetry observed in the WMAP best 
fit power spectrum. Finally we determine the contribution of the modes at very low I to the observed 
asymmetry. 

This paper is arranged as follows. In the next section, a different statistic to understand this even-odd 
power asymmetry is presented. Then, in section 3 we present our results obtained from mock cleaned data 
used to test the effect of foreground residuals. In section [4] we explore the effect of unknown influences on 
the data which might be modulating the primordial signal to induce the observed power asymmetry. In 
section [5] we present our analysis of parity asymmetry in IPSE cleaned temperature data and the cleaned 
maps available using other procedures such as Needlet ILC, etc. In section^ our results from implementing 
different cut at various low— I modes are shown. Finally, we conclude in section [7] 
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2 An Alternate Statistic 



In this section we introduce a different statistic to quantify the parity asymmetry. As we shall see this 
statistic is a more sensitive probe than the one given in Eqs. [5] and [3]. Instead of taking averages of 'V 
even or odd multipoles, we look at local 'I ' power asymmetry. It is defined at0 



where the maximum, l ddi is any odd multipole I > 3 and the summation is over all odd multipoles upto l dd- 
Thus Q{l dd) is a measure of mean deviation of the ratio of power in an even mulitpole to it's succeeding 
odd multipole from one, if it is present in the data. At low— I, since 1(1 + 1)C; ~ constant, statistically we 
expect our statistic to fluctuate about one, like P(l m ax)- 

3 Statistical Significance of parity asymmetry using pure and fore- 
ground cleaned CMBR maps 

We test the WMAP seven year best fit CMB temperature power spectrum^ for anomalous parity asymmetry 
against both pure realizations of CMBR and simulated cleaned maps. The pure CMB sky maps are gen- 
erated as constrained realizations of best fit theoretical CMB power spectrum from AC DM model 2 . The 
synf ast facility of freely available HEALPi^ software (Gorski et al., 2005) was used to produce full sky pure 
CMB realizations at N s id e — 512 of HEALPix's sky pixelization scheme. We then generate five raw maps 
corresponding to each frequency channel in which WMAP makes the observations. We do so by adding the 
pure CMB maps with synchrotron, thermal dust and free-free emission templates from the pre-launch Planck 
Sky Model (PSM)Q (PLANCK Blue Book, 2005) available in each of the WMAP's frequency bands. These 
maps were convolved with appropriate beam transfer functions of K, Ka, Q, V and W bands of WMAP 2 
simulating the raw satellite data. We have also added Gaussian random noise in each pixel using the mean 
rms noise levels in each of the WMAP's frequency channels provided in it's seven year data release (Jarosik 
ct al., 2011). Strong foreground contamination to the observed cosmic CMB signal are assumed to be due 
to galactic synchrotron, dust and free-free emissions. Synchrotron radiation is emitted from relativistic 
electrons in cosmic rays spiraling into the galactic magnetic field. When the dust grains in the interstellar 
medium get heated, they emit radiation due to vibrational mode transitions in infrared frequency which 
is the thermal dust emission. The free-free emission is due to the electron-ion interactions in the ionized 
medium between clusters of galaxies. The full sky simulated raw maps thus generated were cleaned using 
IPSE method as described below. 

3.1 IPSE cleaning procedure 

Here we briefly outline the cleaning procedure we employ to clean the simulated raw maps. The IPSE 
cleaning procedure (Tegmark, de Oliveira-Costa & Hamilton, 2003; Saha, Jain & Souradeep, 2006; Eriksen 
et al., 2007a; Saha et al., 2008; Samal et al., 2010) is a minimum variance optimization method better suited 
for multi-channel CMB observations such as WMAP and PLANCK. It exploits the frequency dependence 
of astrophysical foregrounds received in various detection channels, enabling us to efficiently extract the 

1 While compiling the references to our present work, we learnt that a similar, but not identical, statistic was used by (Land 
and Magueijo, 2005b), which was referenced in the introduction. 
2 Available at http://lambda.gsfc.nasa.gov/ 
3 Available at http://healpix.jpl.nasa.gov/ 
4 http: / /www. planck.fr/heading79. html 
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cosmic signal. The method involves linearly combining various multi-channel maps in multipole space with 
appropriate weights as 

C an = E^#' ( fi ) 

i=l 1 

where af^ an is the clean CMB signal extracted from the raw data, a\ m , acquired from measurements in 
n c frequency channels by linearly combining them with appropriate weights w\. The B\ factors are the 
symmetrized beam transfer functions in multipole space corresponding to an i th frequency channel. The 
weights are computed using the empirical covariance matrix, 

m— — I 

in the formula, 

where eo = (1...1) is a column vector with n c unit elements and W; is also a column vector given by 
(m)j 1 ...m)™ c ) t . The clean power spectrum is then given by, 

/Mclean /n\ 

' ~ TP" 1 ' ( ' 

Further, taking into account the spatial variation of the foreground power across the sky, each map is divided 
into disjoint sky regions and this procedure is applied iteratively in each of these sky partitions. 

Using the IPSE cleaning procedure we generated an an ensemble of 800 cleaned maps. A residual 
foreground bias correction was implemented on each of the cleaned maps by subtracting a bias map estimated 
from these simulations in pixel space (Bennett et al., 2003c). The power spectrum of each of these simulated 
maps was computed using the anaf ast facility of HEALPix and corrected for beam and pixel window effects. 
We used full sky cleaned maps' power spectrum for the range I — [2, 10]. The low— / power can be recovered 
reliably from full sky cleaned maps using the IPSE method (Tegmark, de Oliveira-Costa & Hamilton, 2003). 
For I > 11, we used full sky C; recovered from partial sky map we get after applying a galactic mask 
excluding the heavily contaminated regions in the sky. The KQ85 mask 2 provided by WMAP science team 
in their seven year data release was applied on the simulated maps and we recovered the full-sky C; using 
the MASTER of CMBR or pseudo-C; estimator (Hivon et al., 2002). At low-/, up to / = 32, the WMAP 
best fit power spectrum is estimated from a low resolution ILC map by using the Blackwell-Rao likelihood 
estimator (Larson et al., 2011). For I > 32, they estimated the multipole power using the same pseudo— C; 
estimator that we use to estimate Ci at high—/. These power spectra form the basis of our analysis. 



3.2 Statistical significance 

We next compute the statistical significance using both the statistics, P(l m ax) and Q(l dd), given in Eqs. [3J 
U and Eq. [5] respectively. We used the best fit CMB temperature power spectrum from WMAP's seven 
year data release as reference data. The values of both the statistics for the best fit power are shown in Fig. 
[I]. 

The random chance occurrence probability of getting a P{l m ax) lower than that of the data for l max £ 
[3, 101] is shown in Fig. |2j. For the case of pure maps we used an ensemble of 10,000 simulated CMBR 
maps. We reproduce the results from (Kim &; Naselsky, 2010) using the pure maps ensemble with lowest 
probability at l max — 22. From the 10,000 pure maps we generated, the probability was found to be 0.0013. 
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Figure 1: The even-odd multipole power asymmetry in WMAP's seven year best fit temperature power 
spectrum in the multipole range I — [2, 101] is shown. The asymmetry is computed using both the P{l m ax) 
statistic (lighter curve) and our parity asymmetry statistic, Q(l dd) (darker curve). 

Also, presented in the graph are significances computed using the cleaned maps ensemble. In this case we 
used only 800 simulated maps due to constraints on computational time. As can be seen, the p— values from 
cleaned maps are slightly higher than the probability estimates from pure maps, but are relatively close. The 
foreground power have strong even parity preference and, in (Kim & Naselsky, 2010), the authors speculated 
that the observed asymmetry could be due to over subtraction of foregrounds during foreground reduction. 
But, eventually they ruled out this possibility. Using cleaned maps, we confirm that this is indeed the case. 
Thus, any residual foreground contamination present in the cleaned maps may not induce a particular parity 
preference in the data. The contribution of noise is negligible to the power at low— I. Since, we studied this 
power asymmetry in a wider multipole range, up to Z = 101, and incorporated noise in the simulated raw 
maps, we also conclude that noise cannot cause this asymmetry. From the ensemble of cleaned maps, the 
lowest probability for this parity asymmetry is again found to be at l m ax = 22 with a chance probability 
of 0.13%. This minimum value for significance is beyond 3<r CL and quickly falls below 2a CL by around 
lmax — 40. Beyond l max — 40 it largely stays below 2a. One interesting thing to note is that the P(l ma x) 
curve in Fig. pQ looks wavy, like it was overlayed by some oscillations. It may be indicative of the presence 
of some underlying modulation (see for example Turner, 1983; Martin & Ringeval 2004, 2006; Wang et al., 
2005; Ichiki, Nagata & Yokoyama, 2010). 

With our estimator Q(l dd), we find that, in almost the entire multipole range I = [2, 101], the significance 
of the parity asymmetry lies consistently outside 2a CL in both cases using the pure maps and the cleaned 
maps. These probability estimates are shown in Fig. [3]. The only exceptions are the significances of the 
multipoles 63, 67 and 69 with the cleaned maps ensemble which are marginally inside the 95% CL. Since, 
P{lmax) involves the sum of all even or odd multipole power up to a chosen l max , which is equivalent to 
the mean power up to that l ma x, it appears that their sum is hiding the true significance. As can be seen 
from the plot, we find that for I € [18, 31] the significance of parity asymmetry using Q(l) on WMAP's seven 
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Figure 2: Probability estimates of parity asymmetry seen in the data using the P(l max ) parity statistic, in 
the multipole range I = [2,101]. The significances are computed using 10,000 pure maps and 800 cleaned 
maps, cleaned using IPSE method. We find no significant difference between the two estimates. As can be 
seen, the most significant results occurs at I = 22 for both the cases and is beyond 3er CL. 



year best fit temperature Ci is outside 3cr as estimated from pure maps with minimum at I — 19. In (Kim 
& Naselsky, 2010), l m ax — 22 is specially singled out, for the p— value is lowest at that l max . Here, we see 
that the p— values of Q(l) in the range I = [18, 31] remain close to their minimum. Hence we don't attribute 
any special significance to a particular multipole where Q{1) is minimum, but rather to the whole range 
I = [18,31]. With the cleaned maps simulated set, we find that the p— value curve has slightly rised, but 
only slightly in comparison to that of pure maps. Hence we argue that residuals in the cleaned data cannot 
cause the observed power asymmetry between even and odd multipoles. In the case of cleaned maps, the 
minimum probability for this parity asymmetry is found in the range I = [22,33] using our statistic. 

In our analysis above we used the IPSE cleaning procedure on the simulated maps. It is clearly better 
to use the same procedure for cleaning both the observed data as well as the simulated maps. We do this in 
section 5.1 where we use the ILC cleaning procedure uniformly for the entire analysis. As we shall see the 
results in that case are also consistent with those obtained in the present section. 



4 Unknown foregrounds 

The even-odd multipole power asymmetry we are studying is a point inversion (PI) symmetry violation. So, 
we constructed some templates with explicit PI symmetry breakdown which may induce a power excess in 
odd multipoles and incorporate them in generating our simulated raw maps. In (Dobler & Finkbeiner 2008a), 
an anomalous haze component was found in the WMAP data. It could be that this anomalous haze has 
such asymmetry. Also, there are sub-dominant foregrounds in microwave frequency region which are not well 
characterized yet (Kogut et al., 1996; de Oliveira-Costa et al., 2002; Dobler & Finkbeiner 2008b). However, 
instead of making any such identifications here, we pursue the analysis including these new templates as 
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Figure 3: The p— values of Q{l dd) for WMAP 7 year best fit power spectrum are shown here. As can be seen, 
the data is consistently anomalous by being outside 2a in the whole multipole range I = [2, 101]. For the 
case of cleaned maps, a value less than 1/800 indicates that we obtained no simulated maps whose statistic 
was smaller than that of observed data. 



some hitherto unknown components. A readily conceivable pattern with such a point symmetry violation is 
a hemispherical power asymmetry. This template is shown in Fig. [J at the top. We scale this template 
by a small factor before including them in simulated raw data production pipeline so that it's contribution 
stays sub-dominant and that the simulated raw maps conform visually with the accumulated raw data from 
observations. The results presented in this section were obtained using 350 simulations. Since, we were 
expecting to hnd a lower significance of parity asymmetry in the presence of asymmetric foregrounds, 350 
simulations are sufficient to probe it up to 3cr CL. 

We used this template in two ways, one in which this template is scaled from K-band to W-band of WMAP 
by a frequency dependent power law function and in the other instance as a constant asymmetric foreground 
component, constant in all frequency channels. In the former case, we chose the scaling factor to be 1/15*' 1 
the monopole intensity of synchrotron from PSM at 23GHz (99fiK). It is further scaled to W-band following 
rigid frequency scaling (Bouchet & Gispert, 1999) with a steep spectral index as F{v) = F{v )(vq/v) 2 - & , 
where F{i/q) is the intensity distribution of a foreground component at a reference frequency vq extrapolated 
to another frequency v. We chose a large spectral index so that this effect dies of at higher frequencies (V or 
W bands) where CMB is supposed to be less contaminated by the foregrounds. These templates, generated 
at five frequency bands of WMAP, were added to raw maps and convolved with appropriate beam function. 
In the latter case, where this asymmetric map is added as a fixed power in each pixel across all bands, it is 
scaled by l/30 tft- the synchrotron monopole intensity at 23GHz. The scaling factor is chosen such that this 
excess power will stay lower than the three dominant foregrounds in each channel. The results are presented 
in Fig. [5] using our statistic. 

We find that the significance decreases, but only a little in both the cases. These probability estimates 
are similar to either pure maps or cleaned maps with only the three dominant foreground components from 



7 




Figure 4: Various explicitly power asymmetric templates used in our analysis. These are generated by 
modifying some of the HEALPix routines. The middle and the bottom ones are used in our studies to relate 
North-South ecliptic power asymmetry found by Eriksen et al. (2004a) and the parity asymmetry that we 
are studying here. 
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Figure 5: The p— value plot from the simualted ensemble set generated including our hitherto unknown 
foreground maps with explicit point inversion asymmetry using the Q{l dd) statistic. The significances given 
here are for (1) power law scaled PI asymmetric map, (2) a constant power asymmetric map and (3) a 
foreground component which has explicit power asymmetry in its even and odd multipoles. 



PSM. We also used exponentially scaled prefactors for this template and found no increase in probability, as 
a test for foreground components which may not be following polynomial scaling laws. 

Then, we generated another template which is explicitly asymmetric in power between even and odd 
multipoles. We generated a; m 's with non-zero values for only azo's and zero otherwise and that the asymmetry 
dies (exponentially) with increasing /. The generated map is also shown in Fig. [3], in the middle. Eriksen 
et al. (2004a) found a hemispherical power asymmetry between the north and south ecliptic hemispheres. 
Motivated by the similarity of the that North-South power asymmetry and PI symmetry violation, we explore 
whether there is any relation between these two phenomena. So, we rotate this explicitly even-odd power 
asymmetric map into ecliptic poles. This template is scaled by l/20 th the monopole intensity of synchrotron 
map from PSM and added it as a constant ecliptic dipolar power excess. Here we used a slightly higher 
power to scale this template compared to the earlier constant hemispherical asymmetric template. 

Again we find no change in significance. It only decreases marginally. To find any significant effect with 
this template, we had to add it at unrealistically high levels. With low intensity level, we do not find any 
increase in probability. The significance estimates with our statistic, Q(l), for this dipolar modulation are 
also shown in Fig. [5]- In Fig. [6], we presented the p— value estimates using the P(l m ax) statistic. 



4.1 Multiplicative modulation 

So far we explored the possibility of only "additive" modulations to the data. Now, we also consider 
"multiplicative" modulations (Gordon et al., 2005; Eriksen et al., 2007b; Bunn & Bourdon, 2008; Hanson & 
Lewis, 2009) which can break PI symmetry. The proposed modulation to the CMB is 

AT(h) — 9(n)(l +AX-n), (10) 
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Figure 6: Same as Fig. [5], but for the P(l max ) statistic 

where A is the modulation amplitude and A is the preferred direction. We chose A to lie along the axis of 
ecliptic poles. Thus, we generated a dipole modulation map of Eq. [TU] also shown in Fig. 0], at the bottom. 

Even with an amplitude of A = 0.3 we find that it cannot induce a particular parity preference. The 
results from modulated pure maps are shown in Fig. [Jj. Also shown there are the pure maps estimates for 
comparison. When the same modulation is applied on the raw maps, there is an enhancement in the parity 
asymmetry in the cleaned map, Fig. [7J. This enhancement of power asymmetry in modulated raw maps 
suggest a measurement artifact rather than any thing fundamental to CMB radiation. We point out that we 
are using A = 0.3 which is relatively large amplitude for modulation. Even with such a large amplitude, we 
don't find much change in the statistic's values of modulated pure maps compared to pure maps themselves. 

5 IPSE and others 

We next test the signal of parity asymmetry in maps cleaned by IPSE and several other procedures. In the 
case of IPSE we perform a full sky cleaning of the temperature raw data from WMAP's seven year data 
release. The power spectrum is computed from full sky cleaned map up to I — 10 and used pseudo— C; 
estimator at higher I after applying WMAP's KQ85yr7 mask. Later we estimate the power at low—/ also 
from masked sky. This allows us to determine the how the parity statistic is influenced by masking. The 
Q(l) and P{l m ax) values in the range / = [2,101] for this power spectrum are shown in Fig. [5] and the 
probability estimates are shown in Fig. [S]. These p— values are computed from 10, 000 simulated pure maps. 

We see that the IPSE cleaned data does not show significant power asymmetry. It was surprising to 
find this result, given that the WMAP seven year best fit power spectrum is found to be highly anomalous. 
So, we also applied the two statistics to other cleaned maps available to us. The maps considered here are 
(1) cleaned map from WMAP's five year raw datfl obtained by using the procedure given in Tegmark, de 

6 http:/ /space. mit.edu/home/tegmark/wmap/ or https:/ /www. cfa.harvard.edu/~adeolive/gsm/index. html 
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Figure 7: The p— values obtained after applying a dipolar multiplicative modulation to pure maps and raw 
maps before cleaning is shown here. The parity statistics of modulated pure maps are not different from 
the pure maps themselves. But the modulated raw maps show a slight decrease in the significance of power 
asymmetry. 



Oliveira-Costa & Hamilton, 2003 (2) Needlet ILC map0 of (Delabrouille et al, 2009) using WMAP's five 
year data and (3) Harmonic ILQf| of (Kim, Naselsky & Christensen, 2008), which is also produced from 
WMAP five year data. We note that the maps (1) and (2) are available at resolution of W-band, just like 
IPSE cleaned map. But the Harmonic ILC map is available at 1° resolution. In all these cases the power is 
obtained from the full sky cleaned map up to I — 10 and pseudo— C/ estimator for I > 10, as in the case of 
IPSE. The parity asymmetry statistic values of these maps are also shown in Fig. [8]. As can be seen, all 
these maps give results close to each other, but do not agree with those obtained using the WMAP seven 
year best fit temperature power spectrum. Hence their significances are similar to IPSE cleaned data, as 
shown in Fig. [5]. 

The power spectrum in all the cases analysed in this section is obtained from full sky up to I = 10 and 
masked sky for I > 10. This is in contrast to the WMAP best fit power spectrum which uses masked sky 
over the entire multipole range. Hence it is useful to determine how the results for the maps considered in 
this section change if we use a pseudo— C\ estimator for I < 10 also. In order to estimate parity asymmetry 
using masked sky over the entire multipole range, we applied the WMAP's KQ85yr7 mask on all these 
cleaned maps including our IPSE map. We then obtained the corresponding full sky pseudo-C/ values for 
these maps. The corresponding significance of the parity asymmetry using both the statistics, is shown in 
Fig. [TO]. With our statistic, we find that all the different maps are relatively close to the WMAP's best fit 
power spectrum. The P{l ma x) statistic also shows a higher significance in comparison to the full sky power 
spectra. Hence, we find that, pseudo— Ci estimator, recovered from masked clean maps, reveals the presence 
of anomalous parity asymmetry in these maps. This is most likely due to the fact that the full sky has 
many heavily contaminated regions where the cleaning may not be very efficient. By masking such regions 



7 http: / /www. ape. univ-paris7.fr/APC_CS/Recherche/Adamis/cmb_wmap-en.php 
8 http:/ /www. nbi.dk/~jkim/hilc/ 
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Figure 8: The parity statistic values using Q(l) (top panel) and P(l) (bottom panel) applied to cleaned maps, 
using several different procedures, arc shown here. The maps used are (1) IPSE map, (2) Tegmark's five 
year cleaned map, obtained using the procedure described in Tegmark, de Oliveira-Costa & Hamilton, 2003, 
(3) Needlet ILC five year map and (4) Harmonic ILC five year map. Also plotted are the power asymmetry 
statistic values for WMAP's seven year best fit temperature power spectrum for comparison. All these maps 
show similar levels of parity asymmetry, but do not agree with WMAP's best fit data. 

we hope to get a better estimate of the true power spectrum of the CMB signal. It has been found earlier 
that some of the large angle anomalies disappear in cut sky maps (Bielewicz et al., 2005; Bernui et al., 2007; 
Efstathiou, Ma & Hanson, 2010; Pontzen & Peiris, 2010; Copi et al., 2011). It is therefore encouraging that 
in the present case the signal is enhanced rather than diminished when we use masked sky. 

The fact that a pseudo— Ci estimator gives a higher significance of parity asymmetry in comparison to 
the power spectrum obtained from full sky may lead one to suspect that the process of masking itself might 
generate some signal of parity asymmetry. In order to study this possibility we determine the significance 
of parity asymmetry in WMAP best fit power spectrum by using simulated masked random realizations of 
pure CMB. We use the KQ85yr7 for this purpose. The resulting significance levels for both the statistics 
are shown in Fig. [TT] . For comparison we also show the results for the case when the power spectrum of 
the simulated maps are obtained from full sky. We find that if we use the masked sky pseudo— C; estimator 
for the random samples, the significance level for parity asymmetry is slightly lower for both the statistics. 
Though there is a net rise in p— values due to masking, the relative change is marginal/low and the signal of 
anomalous parity asymmetry is still present. 



5.1 Statistical significance using the ILC procedure for foreground removal 

So far we have used the WMAP best fit power spectrum in our analysis. We computed the statistical 
significance by comparing the statistic for the best fit power with that obtained from the randomly generated 
pure CMB maps as well as simulated, foreground cleaned CMB maps. The simulated CMBR maps were 
cleaned using the IPSE procedure. We have also studied the parity asymmetry in other cleaned maps obtained 
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Figure 9: The p— values for parity asymmetry for the cleaned maps using several different procedures (see 
text). 



using various procedures. In this section we compute the statistical significance using the ILC procedure for 
foreground removal. The ILC procedure is used both for estimating the statistic for the WMAP data as well 
as for cleaning the simulated maps. 

In the ILC procedure, foregrounds are removed by making linear combination of maps at different fre- 
quencies in pixel space. Maps at different frequencies are smoothed to a common resolution of 1° and added 
with suitable weights to minimize the foreground power in the combined map. The details of the ILC proce- 
dure are given in (Bennett ct al., 2003c). In Fig. |12| we compare the power extracted by our implementation 
of the ILC procedure with that obtained by WMAP. We find that the two are in good agreement with one 
another. 

In earlier papers, it has been shown that both the IPSE and ILC procedures are expected to have some 
bias at low The foreground cleaning procedure removes some extra power and hence the extracted signal 
is lower in comparison to the real signal. This effect is dominant at low multipoles I = 2,3. Here we compute 
this bias for ILC using 600 simulations. The extracted bias is also shown in Fig. [T2| . As expected we find 
a negative bias at low— I in power spectrum estimation (Hinshaw et al., 2007; Saha et al., 2008; Chiang, 
Naselsky & Coles, 2009). However we find that the bias is much smaller in comparison to that obtained 
using IPSE. The final power spectrum after removing this negative bias is also shown in Fig. [12] . 

In Fig. [T3] we have shown both the statistics computed for the ILC cleaned map. The corresponding 
statistical significance using the two statistics is shown in Fig. [T3]. The results were presented for the 
WMAP seven year ILC map, the ILC map obtained by us as well as the low— I bias corrected ILC power. 
We find that the statistical significance of all the three maps are comparable to one another. We note 
that the ILC map is reliable on angular scales greater than 10° (Bennett et al., 2003c). It is available at a 
resolution of 1° and so is HILC map of (Kim, Naselsky & Christensen, 2008). So, it will not be meaningful 
to assess the parity preference in that data at high I, even if it shows such an asymmetry, where it's power 
spectrum deviates away from the theoretical CMB power spectrum. 
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Figure 10: The significance using the pseudo-C; estimator. The different maps mentioned in the legend are 
masked using KQ85yr7 mask. 

6 Parity significance with low—/ cuts 

From Fig. pQ, we see that the values of both the statisics, Q(l) and P(l m ax), are much lower at small '/'. 
This is also reflected in Fig. [2] and Fig. [3] where the p— values are found to be relatively insignificant at 
low— I. Also, as mentioned earlier, many studies found that the low— I are associated with various anomalies. 
Hence it is reasonable to assess the parity asymmetry neglecting some of the low— I multipoles. In order 
to get a better insight into this parity asymmetry issue and avoid any "anomalous low— I" concerns, we 
discard some low— I values in this analysis and compute both the statistics for WMAP data and assess it's 
significance. Thus our statistic now becomes, 



Q(lodd) 



la 



l=laut 



^l-l 
Cl 



(11) 



where l cu t is any odd I > 3 and the summation is again over all odd multipoles < l dd- We implement a 
similar /—cut for P(l max ) at low— I in computing P + and P~ . 

The result of applying the two statistics to WMAP best fit power with various low— I cuts is shown in Fig. 
[TS] , As can be seen, both the asymmetry statistics rises closer to one with increasing multipole cuts. With 
different low— I cuts, the p— values at various I in the range / = [l C ut, 101] are computed for both the statistics 
and the results are given in Fig. [TB] and in Fig. [T7]. We see from these figures that the significance of 
parity asymmetry immediately starts decreasing and falls below 2a just by ignoring the first 6 multipoles 
(I = 2, ..,7). This happens with both the statistics. It shows that the dominant contribution to the parity 
asymmetry arises from a few low I multipoles only. It raises the interesting possibility that this might be 
related to the other low-/ anomalies seen in the WMAP data. These might arise from a common physical 
origin. 
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Figure 11: The significance for the WMAP best fit power spectrum using masked sky random realizations 
of pure CMB. The results where the power spectrum of random realizations is estimated from full sky maps 
is shown for comparison. 



7 Conclusions 

We have analysed the signature of parity asymmetry, recently found in WMAP's best fit temperature power 
spectrum, in considerable detail in the multipole range I = [2, 101]. For this purpose, we used the statistic, 
P(l m ax), introduced earlier as well as a new measure, Q(l dd), which appears to be more sensitive. We 
confirm the signal of parity asymmetry at significance level of 3cr. By comparing an ensemble of simulated 
foreground cleaned maps with the WMAP best fit power spectrum, we deduce that the observed parity 
asymmetry cannot be attributed to foreground cleaning or residual foregrounds. Here we use the foreground 
templates from PSM. The PSM may not correctly model some sub-dominant or unknown foregrounds. 
Hence we created some templates which explicitly violate parity and included them in our analysis as both 
additive and multiplicative modulation to CMB. We again find that these cannot explain the observed power 
asymmetry. The level of asymmetry present in data can be obtained by introducing an unrealistically large 
value of these foreground components. Hence we find that we are unable to attribute the observed signal to 
foreground cleaning or residual foregrounds. 

We next tested the presence of this signal of parity asymmetry in several other foreground cleaned maps 
such as the IPSE map, cleaned map using the procedure of Tegmark, de Oliveira-Costa & Hamilton, 2003, 
Needlet ILC map and Harmonic ILC map. We find that these maps also show a signal of parity asymmetry 
provided we use the pseudo— C; estimator after applying a mask in order to eliminate the heavily foreground 
contaminated regions. The significance level for these maps, however, is found to be not as high as that in 
the case of WMAP best fit power. The ILC map also shows parity asymmetry with results closer to that 
obtained with the WMAP best fit power spectrum. 

Finally, we tested the WMAP data for parity asymmetry by eliminating some of the low—/ modes. The 
low— I multipoles are known to show some anomalous results such as low quadrupole power (Bennett et 
al. 2003b), alignment of various multipoles (de Oliveira-Costa et al., 2004) etc. Hence it is possible that 
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Figure 12: Power spectrum of ILC cleaned map from our implementation and from the WMAP seven year 
ILC map is plotted here. The curve below zero of y-axis is the bias in the power spectrum from ILC cleaning 
method. This is computed as an average over 600 simulated pure maps and clean maps generated at 1° 
resolution. Also shown is the best fit theoretical CMB temperature power spectrum for comparison. Our 
implementation and WMAP cleaned ILC map agree with each other. Any difference could be due to our bias 
correction map estimated using PSM. WMAP uses MEM foreground templates for generating the foreground 
bias map. 



the parity asymmetry might also get a large contribution from these multipoles. We found that the parity 
asymmetry disappears by just ignoring the first six multipoles (/ = 2, ..7). Hence we conclude that the low-/ 
multipoles give dominant contribution to the signal of parity asymmetry. It is, therefore, possible that all 
the low-l anomalies, including the parity asymmetry might have a common origin. 
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Figure 16: The p— values of our parity preference estimator, Q(l), as applied to WMAP's best fit power 
with various low—? cuts, using an ensemble of 10,000 pure CMB realizations. The significance falls below 
95% CL just by ignoring the first 6 multipoles. 
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Figure 17: Same as Fig. [TB] but for P(l) statistic, 
disappears by ignoring the first 6 multipoles. 



Here too we find that the significance we found earlier 
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